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Abstract 

The standard map is a paradigmatic one-parameter (noted a) two- 
dimensional conservative map which displays both chaotic and regular 
regions. This map becomes integrable for a = 0. For a / it can 
be numerically shown that the usual, Boltzmann-Gibbs entropy Si(t) = 
~ Si Pi{t) hipi(i) exhibits a linear time evolution whose slope hopefully 
converges, for very fine graining, to the Kolmogorov-Sinai entropy. How- 
ever, for increasingly small values of o, an increasingly large time interval 
emerges, before that stage, for which linearity with t is obtained only 
for the generalized nonextensive entropic form S q (t) = 1 ^jj^ — with 
q — q* ~ 0.3. This anomalous regime corresponds in some sense to a 
power-law (instead of exponential) mixing. This scenario might explain 
why in isolated classical long-range iV-body Hamiltonians, and depending 
on the initial conditions, a metastable state (whose duration diverges with 
1/N — * 0) is observed before it crosses over to the BG regime. 

PACS numbers: 05.20.-y, 05.45.-a, 05.70.Ce 

In his critical remarks about the domain of validity of Boltzmann principle, 
Einstein stressed JlJ that the basis of statistical mechanics lies on dynamics. 
Intensive work has recently been done which is consistent with this standpoint, 
specifically in situations where anomalous effects may arise (j^, 0, 0, B, 0, R]). 
Also, a particularly interesting observation was made in Q, where it was found 
a simple connection between the Kolmogorov-Sinai (KS) entropy rate (the one 
that stems from the properties of mixing of the system) and the statistical 
entropy (the one that originates from a probability distribution) . It was shown 
in fact that, partitionning the phase space in W cells and starting from many 
points within one cell, the time dependence of the usual, Boltzmann-Gibbs 
(BG), statistical entropy 

w 

Si = ~^ j p l \np l , (1) 

i=i 

includes a linear stage whose slope coincides with the KS entropy rate, cal- 
culated, via Pesin equality, using the Lyapunov exponents. Finally, for one- 
dimensional dissipative systems at the edge of chaos (e.g., the logistic map), it 
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was found in |t| that the usual entropy fails to exhibit (with nonvanishing slope) 
such a behavior. Instead, the nonextensive entropy j^] (for a recent review, see 

1 



S « = T — > ( 2 ) 



succeeds for a specific value of q < 1 (q* ~ 0.2445 for the logistic map). We 
remind that the entropic form (||) reduces to (|l|) in the limit q—*\ and that its 
extremum at equiprobability (pi — 1/W, Vi) is given by S q = (W 1 ~ q — 1)/(1— g). 

In Hamiltonian dynamics, one of the most studied models is the standard 
map (see, for instance, PI), also referred to as the kicked-rotator model: 



Xt+i = Vt + TT sin(27rx t ) + x t (mod 1), 

Z7T 



Vt+l = Vt + 7T" sin(27ra; t ) (modi), 

Z7T 



(3) 



where a E M. It corresponds to an integrable system when a = 0, while, for 
large-enough values of a (typically \a\ > 7), the map is strongly chaotic 
It can be easily verified that, for all values of a, this two-dimensional map is 
conservative; it has a pair of (x, y)-dependent Lyapunov exponents, which only 
differ in sign (simplectic structure) and globally vanish when a vanishes. 

In this letter we study the time evolution of both the extensive (q = 1) 
and the nonextensive (q ^ 1) statistical entropies of the standard map (^]), 
focusing our attention on small values of the parameter a, i.e., on those situations 
where the border between the chaotic and the regular regions has a relevant 
influence. For decreasingly small values of \a\ below say 7, one encounters in 
fact an increasingly rich fractal-like structure, characteristic of a large class 
of Hamiltonian systems. Typically, chains of regular islands corresponding to 
elliptic points in resonance condition begin to appear, and between each couple 
of elliptic points of a chain there is an hyperbolic point, thus forming another 
chain of hyperbolic points responsible for the chaotic areas. Around each island 
there is another chain of islands of a higher order and, once again, another chain 
of hyperbolic points between the islands (see, for instance, |Q for details). This 
structure is usually referred to as "islands-around-islands" . If in some sense we 
think at the islands-around-islands structure as the edge-of-chaos region of a 
Hamiltonian system, we verify that for the standard map its relative area of 
influence increases if a approaches 0. 

We know from the theory of chaos that the mixing (herein we use this word as 
a synonym of sensitivity to initial conditions: = lim-A^ro)— >o Aa;(i)/Ax(0), 
for a one-dimensional illustration), is typically linear (£ — 1 oc t) in the case of a 
regular region, and exponential in the case of a strongly chaotic one (£ = Ai£, 
hence £ = e Alt , with the Lyapunov exponent Ai > 0). When a point falls in 
the neighborhood of the islands-around-islands structure, it becomes trapped, 
a phenomenon called "stickiness" (see, for instance, pq]); its motion may be 
thought in this case as a succession of "flights and sticks" . We suggest here, on 
the basis of the results that follow, that the presence of the fractal-like structure 

x We are not interested here in the accelerator-mode islands that appear at various critical 
values of a (see, for instance, j^] and references therein). 
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at the border between regular and strongly chaotic regions produces a power-law 
mixing of the kind (0): 

£ = 4«** = [i + (i-ov*F^ (4) 

(see [[l0| for details on the ^-exponential function), solution of the nonlinear 
equation £ = A 9 «£ 9 , (q* < 1, A 9 « > 0). If we are then to construct a global 
solution representative of the mixing properities of the whole phase space, 
we must implement some averaging process over regions with different kinds of 
mixing. In this case we expect that one or even two crossovers may happen when 
time increases, more precisely between the linear and the power-law mixings, 
and, later on, between the power-law and the exponential mixings. One way to 
model such a behavior is to assume that the differential equation that controls 
£ has the form @ £ = A x £ + (A g . - Ai)f«*. The solution 



£(*) = 



l V | V e (l-g*)Ait 
Ai Ai 



(5) 



presents, in the case < Ai <C A g * , three asymptotical behaviors, namely (i) 

linear (q = regime): £ ~ 1 + A q *t, for < t <C t cross i = nr^lA » » P ower ~ 

i i 9 
law (g = g* regime): £ ~ [(1 - g*)A g »] T ^ T i w 5 ", for t cross i < t < £ cross2 = 

(1 _ g 1 ») Al ; (in) exponential (q = 1 regime): £ - (xJ-)~e Alt , for i > i cross2 . 

With this picture in mind, we explore now numerically the behavior of the 
entropy with time. We start introducing a coarse-graining partition of the map- 
ping phase space by dividing it in W cells of equal size, and we set many copies 
of the system (JVf points) in a far-from-equilibrium situation putting all the M 
points randomly or uniformly distributed inside a single cell. The occupation 
number Mi of each cell i (XL=i = M) provides a probability distribution 
Pi = Mi/M, hence an entropy value. Using then the dynamic equations (||), 
at each step the points spread in the phase space causing the entropy value to 
change. In order to obtain a numerically consistent definition of the probabili- 
ties {pi} when the system spreads at its maximum, we set M > 10 W maxi where 
W m ax < W is the maximum number of cells that the initial data occupy for 
the steps we are considering. Finally, to extract a global quantity on the (x, y) 
phase space, we repeat the calculation setting the cell that contains the initial 
points in different positions chosen randomly all over the whole unit square, 
and we take an average over all the different histories thus obtained. We stress 
here that in taking this average all over the phase space we are consistent with 
a Gibbsian point of view, thus obtaining a statistical description of how the 
system approaches equilibrium. The result of this analysis for fixed a is then a 
single curve of the entropy versus time, for each entropic form S q . 

When the mixing is exponential and we use q = 1, we obtain the behavior 
illustrated in B : a linear increase just before the saturation due to the fmiteness 
of W, with a slope approaching, for W — > oo, the positive Lyapunov exponent. 
If we vary q, we have that for q < 1 (g > 1) the entropy grows with positive 
(negative) concavity before saturation. The exponential mixing is telling us that 
there is an unique value of q, namely q = 1, which produces a linear growth 
in the entropy. We schematize this analysis in Fig. 1. Nevertheless, if some 
phenomenon like the trapping one described in p3| causes the mixing to slow 
down, thus becoming power-law, the same analysis exhibits that the only value 
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of q for which the entropy displays a linear stage with time is q = q* , with 
q* < 1 (sec, for example, [Q). Finally, if the mixing is linear, the specific value 
of q corresponding to a linear stage is q = 0. 

In the case of the standard map the situation is more intricated because 
we have a coexistence of regions with different mixing properties. Fixing, for 
example, a = 5, the phase space presents, macroscopically, two islands and a 
connected chaotic sea. Depending on where we set the initial data, we may face 
a strongly chaotic, or a regular, or even an islands-around-islands region. Our 
averaging procedure does not privilege any of these regions and allows the pre- 
dominant one to emerge. The key point is to perform many different histories so 
that the average stabilizes on a definite curve. The resulting entropy curves are 
then sensible to two different phenomena: the characteristics of the mixing of 
each area (represented by the Lyapunov exponent Ai, or its generalization A 9 * 
in the case of a power-law mixing) and the relative extension of each area with 
respect to the whole phase space (for this reason, in the following we indicate 
the slopes of the linear phases of S q (t) with X q ). If the phase space exhibits 
a predominance of power-law mixing region associated with a generalized Lya- 
punov exponent A g *, with q* < 1, we expect, for time not too large, a linear 
growth for the entropy S q - . Waiting enough time, a crossover to the exponential 
mixing would occur, due to the rapidity of the exponential growth with respect 
to the power-law one, regardeless of how small is the exponential mixing region. 
After the crossover, the linearity with t is obtained for Si (and not for S q *). We 
represent this more complex situation in Fig. 2, in the limit W — » oo. 

Let us present now our computational results for typical values of a. Fig. 
3 represents the case a = 5 and exhibits the predominance of the exponential- 
mixing. After the first step, that here and in what follows we consider as a 
transient, only the BG entropy Si(t) has a linear stage before the effects of 
saturation, while S q has positive concavity for q = 0.7 and negative concavity 
for q = 1.3. For q = 1, our estimation of the slope Ai(5) — 0.96 is different from 
that found in j^] because here we average all over the phase space, including the 
islands. 

We notice that the scale that each entropic form covers, varies widely. To 
perform a linearity analysis between different entropy curves we use the following 
quite sensitive method. First, we fit the relevant piece of each curve (between 
t = ti and t — t2) with a polynomial of the second degree: At 2 + Bt + C. 
Then, we define the linearity coefficient p = A /(At 2 , + Bt^ — At\ — Bt\)\ p 
tests the value of the second derivative of the fitting polynomial, including a 
scale renormalization. Obviously, p = corresponds to linearity, while p > 
and p < correspond to positive and negative curvatures respectively. An 
important remark is that the effect of saturation for a finite value of W decreases 
the curvature of an entropy curve with respect to the case W — > oo. This in 
turns introduces a systematic underestimation of p, which increases when t\ and 
ti increase. 

In Fig. 4 we show how, decreasing a, the first part of the Sy{€) becomes 
concave. After the crossover the curves exhibit once again a linear increase. 
The slope of the linear stage becomes smaller, which suggests that the Lya- 
punov exponent Ai is globally going to 0. In correspondence with the negative 
concavity of the first part of Si(t), So. 3(0 displays a linear increase (Fig. 5), 
with a slope that does not decrease as fast as the one of Si (t) (see the inset of 
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Fig. 4). After the crossover, So.3(*) becomes convex. We notice also that the 
crossover time t cross increases when a decreases (inset of Fig. 5). It is the same 
linearity analysis performed in Fig. 3(d) which leads to q* ~ 0.3. 

Summarizing, we have studied the production of entropy for a well known 
low-dimensional conservative map controlled by the parameter a. It appears 
that, for fixed a ^ and in the limit t — > oo (to be in all cases taken after 
the W — > oo limit), it is the standard, BG entropy Si, which is associated with 
a finite entropy production. However, during the time interval preceeding this 
extreme limit, an increasingly large interval emerges for which it is S q * , with 
q* ~ 0.3, which is associated with a finite entropy production. Consistently, our 
results suggest that, in the limt->oo lifn a ^o ordering, the only regime in fact 
observed is that corresponding to q* ~ 0.3. This fact provides a very appealing 
scenario for (meta) equilibrium thermostatistics in long-range Hamiltonians such 
as that considered in |l5[| . For these systems, a longstanding metastable state 
can exist (preceeding the BG one) whose duration diverges with N — > oo, and 
whose distribution of velocities is not Maxwellian, but rather a power-law. The 
role played by a in our present simple map and that played by 1 /N in such long- 
range-interacting many-body systems, might be very similar, thus providing a 
dynamical basis for nonextensive statistical mechanics B. 
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Figure 1: Schematic representation of S q (t) for different values of q, in the case of ex- 
ponential mixing. Only for q = 1 we have a finite entropy production (limt^oo S q (t)/t). 
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q < q* q = q* q* < 




teross t 



Figure 2: Schematic representation of S q (t) for different values of q, in the limit 
W — > oo and when there is coexistence between power-law (q* < 1) and exponential 
mixings. The thicker lines emphasize the linear regions. 
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Figure 3: Typical S q (t) for a = 5 (a,b,c). 1000 histories were averaged; squares 
correspond to M/10 = W = 708 x 708; circles correspond to M/10 = W = 224 x 224. 
The lines are guides to the eye. (d) Linearity analysis p(q) for the squares, with ti — 2 
and t2 — 6 (see text). The slope of Si(t) between t\ = 2 and t2 — 6 is Ai(5) = 0.96. 





Figure 4: Si(t) for typical values of a (5000 to 7000 histories were averaged); squares 
correspond to M = W = 1000 x 1000 for a = 2, 1.5, 1.1 (M — W — 2236 x 2236 for 
a = 1); circles correspond to M — W — 448 x 448 for o = 2, 1.5, 1.1 (M = W = 
1000 x 1000 for o = l). The lines are guides to the eye. Inset: slopes of Si(t) and 
So.3 (t) (see also Fig. 5) in their linear stage. 
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Figure 5: So. 3 (t) for typical values of a (5000 to 16000 histories were averaged); squares 
correspond to M = W = 2236 x 2236; circles correspond to M = W = 1000 x 1000. 
The lines are guides to the eye. The arrows indicate the crossover time t cross for the 
corresponding values of a, as represented in the inset. 
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